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Abstract 



We study q = 10 and q = 200 state Potts models on dynamical triangulated lattices 
Q^ and demonstrate that these models exhibit continuous phase transitions, contrary to 

the first order transition present on regular lattices. For q = 10 the transition seems 
to be of 2nd order, while it seems to be of 3rd order for q = 200. For q = 200 the 
q phase transition also induces a transition between typical fractal structures of the 

piecewise linear surfaces corresponding to the triangulations. The typical surface 
changes from having a tree-like structure to a fractal structure characterizing pure 
gravity when the temperature drops below the critical temperature. An investigation 
of the alignment of spin clusters shows that they are strongly correlated to the 
underlying fractal structure of the triangulated surfaces. 



1 Introduction 

Important progress was made when it was understood that Euclidean quantum field 
theories can be described as the scaling limit of statistical models on hypercubic lat- 
tices. When a correlation length diverged at a critical point of the statistical model 
one could forget about the underlying lattice structure and define a continuum field 
theory. This point of view added many new aspects to our understanding of quan- 
tum field theory and especially to the concept of renormalization. Two-dimensional 
quantum gravity coupled to matter also allows a lattice regularization, where space- 
time, viewed as a two-dimensional manifold, is represented by a triangulation, i.e. 
approximated by a piecewise linear manifold [3, 2, 4]. Due to the dynamical nature 
of space-time in quantum gravity we cannot confine ourselves to a fixed lattice, but 
have to sum over all possible triangulations of a fixed topology. 1 Here we will only 
consider manifolds with spherical topology. The statistical model is now defined by 
specifying the matter interaction and the observables on an arbitrary triangulation 
and defining the partition function and average value of observables with respect to 
the annealed average over all triangulations. 

As an example consider the g-states Potts model. On a regular two-dimen- 
sional lattice we can define it by 

Z(/3,AO = ]Texp(/3]r^-l] , (1) 

{Qi} \ (m) / 

where 0{ £ {l,...,g}, i denotes a lattice site, {a{\ a spin configuration on the 
lattice and J2(i j) the summation over all neighboring pairs of lattice sites. Finally 
N denotes the volume, which we take to be the number of lattice sites. Let us 
now couple the model to quantum gravity by the use of dynamical triangulations. 
Consider a triangulation T of the sphere and place the Potts spins in the center of 
each triangle. The index i in (1) now labels the triangles. 2 The neighbors of a given 
triangle i are the triangles having a link in common with it. This defines Zj{fi) for 
the given triangulation by formula (1), N being the number of triangles in T. The 
annealed partition function for a finite volume N is then defined by 

Z((3,N)= J2 MP)- (2) 

TeT(N) 

where T denotes a suitable class of triangulations of the sphere. We will later discuss 
in some detail the precise choice of T. 

The g-state Potts models on regular two-dimensional lattices have a critical 
inverse temperature (3 C . Below (3 C the spins are not aligned. Above (3 C the system is 



One could also attempt to sum over all possible topologies, but we will not consider this 
possibility here. 

Of course one could also have placed the spins at the vertices in the triangulation. Historically 
they were placed at the triangles and we follow this convention. 



magnetized in the infinite volume limit. For q = 2, 3 and 4 the phase transition is 
second order and the associated central charge is c = 1/2,4/5 and 1, respectively. 
The q = 2 (Ising) model can also be solved explicitly when coupled to gravity, i.e. 
the partition function (2) can be calculated by use of matrix model techniques. It 
is remarkable that although it is of course impossible to solve the models on some 
fixed random lattice, it is in fact easier to solve the model (2), i.e. the annealed 
average over all random lattices, than to solve (1) on a regular lattice. This might 
be a reflection of the underlying reparametrization invariance of the theory at the 
critical point. One finds that the critical exponents of the models are changed 
relative to the values they have on a fixed regular lattice, and that this change 
agrees with calculations performed directly in a continuum framework, considering 
two-dimensional quantum gravity coupled to conformal field theories with central 
charges c < 1. 

It is still something of a mystery what happens when conformal field theories 
with c > 1 are coupled to quantum gravity. The continuum approach leads to 
complex critical exponents and it is presently not known how to repair this. Recently 
progress has been made in understanding the physics for c > 1 using the lattice 
approach [5, 6, 25, 27, 17]. Especially the limit c — > oo is now reasonably well 
understood. In addition it can be proven that the low temperature limit of the 
q — > oo Potts model agrees with the low temperature limit of the multiple Ising 
model when the number of copies n of the Ising models goes to infinity [26]. Since 
the central charge of the multiple Ising system will be ra/2 we see that the large q 
limit of the Potts models provides information about the large c limit. 

At first sight this might be surprising since it is well known that the g-state 
Potts models has first order transitions for q > 4 on regular lattices while the 
multiple Ising models on regular lattices still have second order transitions. Second 
(or higher) order transitions allow us to think about divergent correlation lengths 
and associated continuum field theories. In the case of first order transitions we 
have no continuum field theory we can associate with the spin model and although 
a statistical model like (2) is perfectly well defined, we have a priori no continuum 
field theory for which we can say: (2) presents a regularized version of this field 
theory coupled to 2d quantum gravity. However, the study of critical properties 
of discretized 2d quantum gravity shows many examples of this kind. We can just 
mention the so-called multicritical matrix models. One observed that they gave rise 
to new critical behavior of the statistical system, but only later the corresponding 
conformal field theories which led to the critical exponents after coupling to gravity 
were identified as the (2, 2m- 1) minimal conformal field theories. From this point of 
view it seems justified to take the following attitude: if the annealed average of a q- 
state Potts model shows non-trivial critical behavior, i.e. exhibits a continuous phase 
transition rather than a first order transition, the corresponding critical behavior 
might indeed reflect the behavior of a conformal field theory coupled to gravity. 
The weak coupling results for q — > oo seem to support this idea and they indicate 



in addition that the behavior is associated with conformal field theories with c > 1 
coupled to gravity for q sufficiently large. This should serve as a strong motivation 
for a detailed study of the properties of the g-state Potts models coupled to 2d 
gravity. In addition the models have course some interest as statistical systems with 
new non-trivial critical behavior. 

In this paper we will show that the critical behavior of the q = 10 states 
Potts model and the q = 200 states Potts model indeed changes from being trivial 
(i.e. first order) to non-trivial, namely second order for q = 10 and third order for 
q = 200, when we change from a regular lattice to the annealed average (2). 

The rest of the article is organized as follows. Section 2 summarizes some 
recent predictions concerning large q and c models. Section 3 motivates the use of a 
particular class of triangulations T. In section 4 we determine the order of the phase 
transition for the two Potts models coupled to gravity. In section 5 we discuss the 
back-reaction of the g-state Potts models on the geometry, defined by the piecewise 
linear manifolds corresponding to the triangulations, and we discuss the correlation 
between geometry and matter, more specifically between so-called baby universes 
and spin clusters. Finally section 6 summarizes our present understanding of the 
g-state Potts models coupled to 2d quantum gravity. 

2 Expectations 

Recently the first detailed predictions on the behavior of two-dimensional quantum 
gravity in the "strong coupling" regime have become available. In particular, the 
oo-state Potts model and a class of c = oo models (including, for example, an infinite 
number of Ising species on the surface) have been solved exactly on a random surface 
using matrix model techniques. The biggest surprise was that all of these models are, 
up to some trivial transformations of the temperature, identical — and this inspite 
of the well-known fact that on a fixed two-dimensional lattice the Potts model has 
first-order, rather than continuous, phase transitions for q > 4. 3 The reason for this 
unexpected equivalence is that, roughly speaking, the q = c = oo model is a mean 
field theory: each state, each point in the (discrete) target space (the complete graph 
on q points for the Potts model, an ra-dimensional hypercube for n Ising species) 
has so many neighbors that it feels only them, rather than the global structure of 
the target space. 

The qualitative picture at q = c = oo is as follows. At high (3 there is the 
usual ordered, pure gravity phase with j s = —1/2. In the low q, c models this phase 
is bordered by the critical point (third-order) where j s jumps to some higher value, 
followed by the low-/3 phase — another phase, this time disordered, of pure gravity. 
In the case of q = c = oo, however, at low (3 the surfaces degenerate into trees, the 
nodes of which are spin clusters; in this phase the clusters are small and the trees 



This explains the relative neglect of the random-surface Potts model at high q. 



are large, and therefore we have what is commonly called branched polymers, with 
j s = +1/2. At the multicritical point between the two phases we have j s = +1/3, 
and the transition is third-order, with a = — l. 4 Thus for the oo-state Potts model 
the transition is softened from first to third order by the coupling to gravity, while 
for the c = oo models the order is promoted from second to third. 

A direct approach to the branched polymer regime, one that avoids use of 
matrix models and deals directly with the triangulated surfaces, is possible [5, 6, 
17]. Starting with theoretical arguments for branched structure at large c, as well 
as computer observations [9], one can postulate such a structure as a kind of a 
"polymerization" of an ordinary random surface model by restricting the boundaries 
of spin clusters to have minimal length. In this way, by assuming the spin cluster 
geometry of the q = c = oo surfaces, one arrives at the same picture of their 
critical behavior as via the exact matrix model solution, with the same values of the 
exponents j s and a. Finally, a novel combinatorial method [19] proves that trees 
are dominant in the c — > oo limit, and finds j s = +1/2 in the low-/? regime. 

What about the Potts model with finite ql We know that at small q, the 
low-/? regime should be a disordered, pure gravity phase rather than a branched 
polymer. One way to get there is to expand the q = oo matrix model in inverse 
powers of q [27]. The first non-trivial correction gives a new structure: at the vertices 
of the trees are no longer spin clusters, but complex graphs of spin clusters with 
(arbitrarily many) loops, structures which we will call "galaxies." By adjusting the 
coupling constants one can now vary three things: the area of the spin clusters, the 
size of the galaxies (i.e., the number of clusters inside the galaxies), and the size of 
the trees (the number of galaxies inside the tree). 5 

In this approximation, for large enough q the critical behavior is as follows: 
there is the ordered, pure gravity phase (j s = —1/2) at large (3, a phase of large spin 
clusters 6 ; there is still a third-order transition (a = —1) at some /3 C 2 to a branched 
polymer phase of large trees that has j s = +1/2; but now, for any finite q at (3 = /3 c i 
there is another third-order transition (a = —1) to a phase where the galaxies are 
large and j s = —1/2 (once again), and which extends all the way down to (3 = 0; 
finally, at both transitions j s = +1/3. This new large galaxy phase can be identified 
as an disordered phase — precisely what is missing at q = oo. As q decreases from 
infinity, /3 c i and /3 C 2 approach each other: the branched polymer, mean field phase 
gets eaten away, precisely as expected; until some q c at which the branched polymer 
phase disappears altogether, leaving us with just the picture we know to be true at 
low q: two pure gravity phases, an ordered and a disordered one, separated by a 
phase transition. Fig. 1 is a schematic two-dimensional slice of this phase diagram. 7 



It should be noted that here, as well as in the Ising model, there are no logarithms in the free 
energy and — as a is an integer — no poles in any of its derivatives: there is simply a discontinuity 
in the third and all higher derivatives at the critical point. 

Subsequent refinements to this approximation should not add any new essential features. 

A "large" quantity in this context means that some of its moments diverge. 

The third dimension that has been suppressed is the cosmological constant: in simulations of 
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Figure 1: Schematic phase diagram of possible phase diagram of the Potts model at 
large q. (I) is the ordered (pure gravity) region of large clusters, (II) is the branched 
polymer region of large trees, and (III) is the disordered (pure gravity) region of 
large galaxies. 



The exact critical values of the parameters of course depend on the truncation 
used in the above approximation (see [27] for details). Two methods were used in 
[27] that gave q c ~ 120 ± 20. The same two methods give /3 C 2 — fi c \ ~ 0.1 ± 0.05. 
The asymptotic behavior of /3 C 2 is ^ log A + 0(A~ 1 ) 1 where A = q— 1. (3 c i is harder 
to calculate, but it seems that for A large enough (3 c i ~ |logA. 



3 The use of tadpoles 

It is by now well established that Monte Carlo simulations may serve as an effective 
tool for extracting the critical behavior of systems like (2). Let us here only mention 
that the Monte Carlo simulations involve two types of updating: we have to update 
the fluctuating random triangulation and we have, for a given random lattice, to 
update the matter fields. The updating of the triangulation uses the so-called flip, 
where two adjacent triangles, viewed as a square with the common link as the 
diagonal, is replaced by the same square but made out of the two triangles formed 
by using the other diagonal in the square (the diagonal is "flipped"). For a given 
triangulation the updating of the spins is performed by the Swendsen-Wang cluster 
algorithm. In most simulations used until now the class of allowed triangulations has 
been characterized as combinatorial two-manifolds. The two requirements satisfied 



large surfaces, this can of course be considered tuned to its critical value. 



by such triangulations are that any two vertices have at most one link in common 
and that all vertices, which are not on the boundary, belong to at least three links. 
It can be convenient to relax this concept of a combinatorial manifold and consider 
all triangulations which can be obtained by successive gluing of triangles along their 
sides to form a closed surface of spherical topology, i.e. so that Nt — Nt + Ny = 
X = 2, where Nt, Nl and Ny denote the number of triangles, links and vertices in 
the closed surface after the gluing. For such triangulations we can have one-loops 
(i.e. links which start and end at the same vertex) and two-loops (i.e. two vertices 
joined by more than one link as mentioned above). To each triangulation we can 
associate a dual (f> 3 graph; the difference between the (f> 3 graphs which correspond to 
combinatorial two-manifolds and the <f) 3 graphs which correspond to the spherical 
triangulations obtained by the unrestricted gluing is that the last set of graphs 
contains tadpole and self-energy diagrams which are associated with the one- and 
two-loops. 

Whether or not one- and two-loops are included should not effect the scaling 
limit of the theory, and we know that it does not in the cases where the theories 
can be solved explicitly, i.e. for central charge c < 1. The non-universal features 
are affected, though. Especially the magnitude of finite-size effects can be different. 
This is important from a practical point of view and we have observed this difference 
both in strong coupling expansions and Monte Carlo simulations. In both cases the 
inclusion of tadpoles and self-energy diagrams significantly reduces the finite-size 
effects. 

Let us first discuss the strong coupling series. We can expand the free energy 
up to order 14 in the cosmological constant g, and use a [2,1] Pade approximant 
on the ratios of successive coefficients, which should behave as g~ 2 (n/n — 2) 7s_3 . 
For pure gravity, where j s = —1/2 for spherical graphs, we find j s = —0.503 when 
we include all graphs, but j s = —0.516 when we exclude tadpoles and self-energy 
diagrams. The relative advantage of tadpoles and self-energy diagrams becomes 
more marked if, say, an Ising model (i.e. a q = 2 Potts model) is coupled to the 
surface. In fig. 2 we plot j s as a function of c = e~ 2 ^ using the Pade approximation: 
the benefit of tadpoles and self-energies is evident. Other Pade approximants and 
other extrapolations like the ratio method give similar results. 

In order to test this effect directly on much larger surfaces we have performed 
Monte Carlo simulations with high statistics (10 6 sweeps) on the Ising model at the 
critical values (3 C for both classes of triangulations . The results are shown in fig. 3. 
Finite-size scaling is used to extract the critical indices related to the magnetisation 
M and the susceptibility \- The comparison is made easier by the fact that the 
exact critical coupling for the Ising model is known for both types of triangulations, 
/3 C = i log(108/23) if one and two-loops are excluded and fi c = \ log((2\/2 — l)/27) 
if they are allowed [15, 16]. 

The simulations are performed on triangulations ranging from N = 20 to 
N = 20480, and from the magnetisation M and the susceptibility \ finite size scaling 
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Figure 2: j s as a function of temperature for the Ising model, using strong-coupling 
expansion. Solid curve is with tadpoles and self-energies included, dashed curve is 
with them excluded. 
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Figure 3: The extracted values of (a) (3/vdu and (b) y/isdu (measured at fi c ) 
versus lattice size for one Ising model coupled to gravity. The data points with 
circles are obtained for triangulations which include one- and two-loops and the 
data points with squares are obtained for triangulations where the one- and two- 
loops are excluded. The solid lines represent the analytical values (3/vdu = 1/6 and 
y/vd H = 2/3. 



allows us to extract the combinations (3/vdu and y/isdu of the critical exponents. 
(3 and \ are the critical exponents of the spin and the susceptibility while v is the 
exponent for the correlation length, dfj denotes the Hausdorff dimension of the 
ensemble of piecewise linear surfaces defined by the triangulations. 

From the finite size scaling relations it follows that (3/vdu and y/isdu in 
principle can be determined from two consecutive lattice sizes N and 2N, i.e. 

Z 3 (iY) = logM(2A0-logM(AQ 7 (y) = logx(2AQ - logx(N) 



vdu log 2 vdu log 2 

In fig. 3 we plot the extracted values versus lattice size. It clearly demonstrates how 
the finite size corrections are diminished if loops of length one and two are included. 
Another motivation for using the triangulations with one- and two-loops 
comes from the conjecture [7] that the transition to the magnetized phase takes 
place only after the surfaces have developed a pronounced tree-like structure where 
the individual branches are already magnetized before the transition, while the total 
magnetisation is still zero. The transition to the magnetized phase is only the final 
alignment of the spins of the individual baby universes. Allowing one- and two-loops 
makes it easier to detect how the location of spin clusters fits to the baby universe 
distribution as the fractal structure becomes more pronounced. 

4 Order of the phase transitions for q = 10 and 200 

4.1 The simulations 

We now turn our attention to the simulations of the g-state Potts models for q = 10 
and q = 200. As for the Ising model we use the flip algorithm to update the 
triangulations and the Swendsen-Wang cluster algorithm for the spins. The lattice 
sizes were N = 250, 500, 1000, 2000, 4000 and 8000. The main problem in the 
simulations is that the autocorrelations seem to grow linearly with q, presumably 
due to the increased branching of the surfaces. 8 This forced us to perform very 
long runs in order to get good statistics. In the vicinity of the phase transition we 
performed approximately 5 X 10 6 and 2 X 10 7 sweeps for each fi value for the q = 10 
and q = 200 state Potts model, respectively. With the autocorrelations present this 
resulted in between 150.000 and 400.000 independent measurements for q = 10, and 
40.000 to 60.000 independent measurements for q = 200, depending on lattice size. 

4.2 Cumulants 

It can be difficult to distinguish numerically between a continuous and a first order 
phase transition. Many methods have been proposed [14] which either depend on 



It would be very interesting to use the new algorithms developed [21] to deal with surfaces of 
pronounced fractal structure on this problem. However, the algorithms were not available when 
this project was started. 



the presence of hysteresis at 1st order transitions or use finite size scaling. All of 
the methods have their drawbacks. Hysteresis can be hard to detect in simulations 
and if one drives the system too fast pseudo-hysteresis can also be seen in second 
order transitions and one quickly enters into discussions of weak 1st order transi- 
tions versus strong 2nd order transition. Finite size scaling only applies when the 
system size is larger than the characteristic correlation length at the pseudo-critical 
transition point and is based on a number of assumptions which are not necessarily 
satisfied. Especially in the case where we consider matter systems coupled to quan- 
tum gravity one should consider the assumptions with some suspicion: what is the 
system size when we integrate over fluctuating geometries of fixed volume? Recall 
the basic equation of finite size scaling: 

L d f L (t) = F(t»L), t = \(3-(3 c \ (4) 

where L is the size of a regular lattice, / the bulk free energy density, (3 the inverse 
temperature, d the dimension of the the lattice and v the critical exponent of the 
correlation length £(/?) ~ l/|/3 — (3 C \ U In our case we have only N, the total volume. 
We should presumably make an identification N = L dH where djj is the Hausdorff 
dimension of our ensemble of manifolds, and we should presumably View £(/?) as some 
average geodesic distance on the same ensemble, but it is clear that the foundation 
for finite size scaling is not as firm as on regular lattices. We have nevertheless chosen 
to use finite size scaling for two reasons: as was clear from the data presented in 
the last section it seems to work very well for matter fields with central charge 
c < 1 coupled to gravity so it is not a vain hope that it will work well for other 
models. Secondly it is the method which allows us to extract critical exponents in 
the cleanest way. 

Let us now look at scaling of the specific heat Cy and the energy cumulant 
V L 

_<£ 4 > 

(E 2 

In the case of a first order phase transition the specific heat will diverge and in 
finite systems its maximum will grow linearly with lattice size, Cv(max) ~ N. In 
fig. 4a and b are shown the measurements of the specific heat both for q = 10 and 
q = 200. The maximum increases very slowly with lattice size, much slower than is 
expected for a 1st order transition. These statements are quantified in fig. ha where 
Cy(max)/N is plotted versus TV -1 for both models. If we were dealing with first 
order transitions this quantity should approach a non-zero constant, but instead it 
seems to go to zero. 

The energy cumulant Vl provides even more conclusive evidence. For a 
continuous phase transition it vanishes in the infinite volume limit for all values of 
(3, while at a 1st order transition fluctuations between coexisting phases prevent it 
from becoming zero in the critical point [12, 13, 18, 20]. In that case Vl reaches 
a non-zero minimum value which scales with N. In fig. Ac and d we show Vl 
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Figure 4: The specific heat Cy, for various lattice sizes, for (a) q = 10 and (b) 
g = 200 state Potts models coupled to gravity. In (c) and (d) the energy cumulant 
Vl for the same models is shown. 
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Figure 5: (a) Cv(max)/N versus 1/N for q = 10 and 200. Linear fits show that 
Cv(max)/N — > as N — > oo. (b) The maximal slope of Binder's cumulant versus 
lattice size. The linear fits yield vdfj- (c) and (d): the pseudo critical points (3 C (N) 
versus N -1 '" H , both obtained from Cy(max) and from intersection of BC. 

for different lattice sizes for q = 10 and 200. In both cases it approaches zero as 
N increases and no evidence of the characteristic behavior of a 1st order phase 
transition (a persistent dip) is seen. 

Taken together these results provide strong evidence in favor of continuous 
transitions. As already emphasized there are no compelling arguments why a 1st 
order phase transition on a regular lattice stays first order after we perform the 
annealed average (2) over random triangulations. On the contrary, for the Ising 
model it is known that the phase transition is weakened from 2nd to 3rd order, and 
numerical simulations show that the same is the case for the 3-state Potts model. 
It is true that in earlier simulations on the 10-state Potts model coupled to gravity 
it was claimed that the transition was 1st order [11]. But that was concluded from 
the qualitative behavior of Binder's cumulant, which is not a very good indicator to 
the order of phase transitions. In fact our measurements are in excellent agreement 
with those in [11], but we must draw a different conclusion. 
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4.3 The critical temperature 

The next step is to locate the critical temperature. It can be done in two ways by 
finite size scaling. One can either use the location of the peak in the specific heat 
(3 C , which scales as 

I&-/3J- A" 1 ^, (6) 

or look at the intersection of Binder's cumulant for different lattice sizes. 

Let us first consider the specific specific heat. As we do not have many lattice 
sizes at our disposal it would be impossible to make a three-parameter fit to (6). 
Fortunately it is possible to determine vdn directly, as the maximum of slope of 
Binder's cumulant scales as 



(dBC_\ 
where Binder's cumulant is defined as 



'jdp 



™{lM-)'~ NH ' < 7 » 



BC - {Mi 



(M 2 ) 2 v ' 

In fig. 56 we show the lhs of (7) on a log-log plot and the corresponding linear fits. 
The values of vdn obtained in this way are listed in table 1. Knowing vdn we can 
now make a linear fit to (6) and obtain (3 C (/3 C (CV) in table 1). 

Next we turn to the use of Binder's cumulant (8) to locate the critical point 
of the two multiple Potts models. As is well known the curves of BC intersect at 
(3 C up to finite size corrections. Our statistics are so good that we can actually see 
how the intersection point moves with increasing lattice size. It is to be expected 
that the intersection point moves towards the true critical point in the same way as 
the maximum in the specific heat, i.e. as given by (6). This is shown, together with 
(3 c (Cv(max)), in fig. 5c and d. Extrapolating to A = oo yields the critical points 
shown in table 1. The two methods for determining (3 C are in fairly good agreement, 
the slight discrepancy in the values is presumably being due to the fact that larger 
lattices are needed before the methods converge completely. 

It should finally be mentioned that the critical temperature we find agrees 
fairly well with the asymptotic prediction [25] fi C2 = \ log q + 0(q~ 1 ) x 2.65. 

4.4 Critical exponents 

To determine the magnetic critical exponents (3, j and a we apply finite size scaling 
at the critical point (taking (3 C as the weighted average of /3 C (CV) and (3 C (BC)). From 
the scaling of the magnetization M and its derivative with respect to (3 (which scales 
as A^ 1- ^' v H ) we get estimates for fi and vdjj . ^From the scaling of the magnetic 
susceptibility we get y/isdij. Finally, and most interesting, the specific heat allows 
us to estimate a, both from the scaling of its maximal value and at the critical point. 
All these exponents are shown in table 1. The fits to the specific heat are shown 
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2.52(12) 


1.412(1) 
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+0.02(4) 
-0.83(6) 


+0.02(3) 
-0.80(5) 



Table 1: The extracted critical points and exponents for the q = 10 and 200 state 
Potts models coupled to gravity. The values are obtained using finite size scaling. 
For a we have values both from the scaling of (a) Cv(max) and (b) CV(/3 C ), and 
for vdfj from the scaling of (a) max{dBC/d(3} and (b) M and dM/dfi. The main 
contribution to the errors comes from the location of (3 C . 

in fig. 6. For q = 10 a is very close to zero and we have thus tested whether the 
specific heat fits better to a logarithmic divergence 



Cv(N) = c 1 + c 2 log(N). 



(9) 



This fit is also shown in fig. 6 and turns out to be equally good as the one with 
a = 0. 

4.5 Comments 

What is the interpretation of these results? The first thing to remember is that for 
the q = 4 state Potts model coupled to gravity we have (3 = 1/2, j = 1, a = and 
udfj = 2, i.e. consistent with our measurements for the 10-state Potts model. It is 
natural to conjecture from this coincidence that we have an interval of g's starting 
at q = 4 and extending beyond q = 10 where we have identical critical behavior 
after coupling to gravity. We will discuss this remarkable possibility later. 

At some large value of q the the phase transition changes character, a has 
a value close to —1 for q = 200. This value agrees both with [27] and [6], indicating 
that we might have reached the large q limit in which the magnetic transition takes 
place between surfaces degenerated to branched polymers (in the phase with no 
magnetization) and pure gravity (the low temperature phase where the system is 
magnetized). We will later show that the entropy exponent for the surfaces gives 
some support to this picture. 

Whether this change of critical behavior happens at some critical value of 
q, or whether there is a gradual change from one behavior to another cannot be 
established from our data. 
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Figure 6: The measured values of Cy, both the maximal value and in the critical 
point (3 C , for (a) q = 10 and (b) q = 200 state Potts model. Included are fits to 
Cy = c + ciN a l vdH and, for q = 10, an equally good fit to Cy = c + c\ log N 

5 The fractal structure of surfaces and spin-clusters 

5.1 Expectations 

Let the partition function of the annealed average of the Potts models over random 
triangulations be given by (2). If we perform the summation over spin configurations 
we can write the partition function as 

Nls (fJ)-3 e i, c (fJ)N 



Z(0,N) 



N 



oo. 



(10) 



This formula defines the string susceptibility exponent J S (P). It is tacitly assumed 
that the partition function can be written in this form and it is well known that it 
is the case for q = 2 and q = 3 where the models can be solved analytically. 

Recall the scenario of the critical behavior of the large q Potts model pre- 
sented in sec. 2. We have tried to substantiate this picture in the case of q = 200. 
There are three questions to address: 

1. Can we observe the branched polymer phase for some temperature range just 
above the transition from a magnetized phase to a phase with total magneti- 
zation equal zero? 

2. Can we observe a second phase transition at higher temperature from the 
branched polymer phase and back to some pure gravity phase? 

3. Can we observe the magnetisation of baby universes and the correlation be- 
tween spin clusters and baby universes? 

5.2 Search for branched polymers 

Let us first look for the branched polymer phase. The exponent j s is a good signal. 
It should be equal 1/2 if the surfaces are degenerate into genuine branched poly- 
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Figure 7: Measurements of j s and r ma:E for g = 200 state Potts model coupled to 
gravity. 

mers. Another quantity, suggested in [7] is the maximal "link" radius (r max ). More 
precisely we define a distance dij between two vertices i and j in a triangulation 
as the length of the shortest path along links between the vertices. r max (T) is now 
defined by 

Ev" T maxj{rffj} 



iax\^- ) 



V T 
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where Vj denotes the number of vertices in the triangulation T. (r max ) is now the 
average value of r max (T) in the ensemble of triangulations defined by the partition 
function Z(N,(3) from (2). Although (r max ) is not a universal quantity like the 
susceptibility exponent j s it reflects very well the geometry of the surfaces and it 
is known that in the case of multiple q = 2, 3 and 4 state Potts models coupled to 
gravity it shows a pronounced dip on the high temperature side of the spin transition, 
the depth increasing with the total central charge of the system. 

In fig. 7 we plot these quantities for q = 200 and a wide range of (3. The 
values of j s are determined from the counting of the so-called baby universes. For 
an ensemble of surfaces given by the partition function (10) the distribution of baby 
universes with area B < N/2 will be given by 9 



(12) 



n{B) ~ NB ls ~ z ll-B/N) 



j s is now determined by fitting the measured distribution to (12). The measurements 
of j s may indicate a branched polymer phase starting at (3 ~ 2.4 and extending to 



As we allow triangulations with one-loops and two-loops we can just as well extract -/ s from 
the distributions of baby universes with neck length one and two as well as length three (as was 
used in the original measurements [8]). We performed simulations on pure gravity in order to test 
this. The results for surfaces of size N = 1000 were: y s (l = 1) = -0.501(6), y s (l = 2) = -0.502(8), 
j s (l = 3) = —0.498(22). We conclude that there is universality for n(B, I) as long as the loop length 
I is very small. 
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the magnetisation transition. We have a rather broad peak value around j s ~ 0.4. 
This broad peak is in sharp contrast to the very narrow peaks one observes for spin 
systems which at the critical point corresponds to a conformal field theory with c < 1 
(i.e. here q = 2 and 3) and it gives some indication that there might be a plateau 
where j s / —1/2. (For q = 2 and 3 it is known that j s is only different from the 
pure gravity value —1/2 precisely at the transition point). We do not get the precise 
value j s = 1/2 for branched polymers. It is not clear to us whether this is finite 
size effects or simply the fact that we are below the critical value q c above which we 
should enter the branched polymer phase. There is the following theoretical reason 
for believing that we observe finite size effects: from general arguments spin systems 
should either have j s < or j s = l/(ra+ 1) [17], and probably N = 1000 is not very 
large for q = 200. Unfortunately it is very time consuming to enlarge the system 
for q = 200, as already mentioned. 

When we look at (r max ) it is important to keep in mind that it is not a 
universal quantity. It might change quite a lot without reflecting a corresponding 
change in the critical properties of the system. In fact it decreases continuously 
from the high temperature limit, where it agrees with the value for pure gravity to 
a minimum at (3 ~ 1.75. After this it increases rapidly to its maximum close to (3 C , 
corresponding to the magnetisation transition. It does not fit in a very convincing 
way with the behavior of j s , but one could argue that it is not inconsistent with the 
following picture: at infinite temperature the spin system has a life independent of 
the underlying lattices, which therefore have the geometry of pure gravity. As the 
temperature decreases the interaction between lattices and the spin system increases 
and the lattices are deformed towards branched polymers 10 . Only deep into this 
phase j s will then jump from —1/2 to 1/2. This suggests that one should look for 
the second transition suggested in [27] in the region 1.75 < /32.5. 

5.3 Search for the second transition 

According to the conjecture in [27] the second transition, where j s jumps back from 
1/2 to —1/2, while taking the value —1/3 at the transition point should be a 3rd 
order transition. From the point of view of the fractal geometry of the surfaces it 
is a kind of mirror transition of the magnetisation transition, but taking place at 
somewhat higher temperature. Clearly it cannot be a spin transition in the usual 
sense, so the only general signature we have is that the critical exponent of the 
specific heat a should be —1, i.e. there should be a discontinuity in the derivative 
of the specific heat. We have scanned values of (3, especially the region suggested 
above, but have found no trace of a transition. 



More precisely the dual lattice will acquire a conventional branched polymer structure since 
the spins are placed on the dual lattice. This is precisely what gives a decreasing value of (r„ 



' max / 



for the triangulated surfaces [7]. 
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5.4 Spin clusters 

Let us now turn to the correspondence between spin clusters and the fractal struc- 
ture. From [25, 6, 27] the magnetisation transition for large q or large c arises by 
alignment of already magnetized baby universes. To which extent is that true? To 
answer that question we propose two new observables: 

• The ratio between number of boundary links 6; (links separating different spin 
clusters) and number of spin clusters Nc- 

W=*±±. ,13) 

• The number _3 of sets of three clusters that all share boundaries with each 
other which we shall call triplets (normalized by the number of clusters). 

Let us explain why these observables tell us something about the alignment 
of spin clusters. First W. It is always larger or equal than one and, more important, 
it can only be equal to one in two situations: (i) if the system is magnetized and 
hence there is only one cluster (no boundary links), (ii) if all the clusters are glued 
together in a treelike structure separated by boundaries of length one. Similarly for 
A3: It can only be zero in the two situations described above. It is however not as 
sensitive to the neck size separating the baby universes. So if the baby universes are 
completely magnetized before the phase transition we should observe W ~ 1 and 
A 3 ^0. 

Fig. 8 shows the measurement of these quantities for q = 4, 10, 50 and 200 
state Potts models coupled to gravity and, for comparison, for the same models on 
a regular square lattice. 

Let us start with a discussion of W (fig. 8 (a))- On a triangulated surface each 
spin has three neighbors so that for very high temperatures, where the spins fluctuate 
independently, each spin will form its own cluster (for q large enough) and every link 
is a boundary link. The ratio between links and spins on a triangulation is 3/2, which 
is what we measure in the (3 — > limit. For (i ~^> (i c there is only one cluster, hence 
W = 1. What is interesting is that the value around the magnetization transition 
decreases as q increases. It starts at W ~ 1.5 for q = 4 and reaches W ~ 1.05 
for q = 200. These values do not change notably with lattice size. To distinguish 
between situations (i) and (ii) we measure the number of spin clusters around (3 C . 
For q = 200 N c ~ -/V/4, which shows that W ~ 1 because the baby universes are 
magnetized individually (in different direction) while the total magnetization is still 
zero. In other words, the graph of spin clusters is large and increasingly tree-like at 
the critical point as q — > 00. 

Very different behavior is seen on a regular lattice. As before the spins fluctu- 
ate independently at high temperature and we have W ~ 2 (the ratio between links 
and spins on a square lattice). For high (3 we get W = 1 due to total magnetisation. 
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W 1.6 




Figure 8: Measurements of W and A3 for q = 4, 10, 50 and 200 state Potts models 
coupled to gravity ((a) and (c)) and, for comparison, on a regular lattice ((b) and 
(d)). The values are plotted versus the reduced temperature r = (/3 — (3 C )/(3 C . 
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But close to the phase transition the value of W is very large as there are large and 
complicated clusters (and which cannot decrease the boundaries by deforming the 
lattice). Furthermore the sharp jumps in W at (3 C show the strong 1st order nature 
of the transitions for q > 4. 

Exactly the same picture is seen for the number of triplets A3. For randomly 
triangulated surfaces the value of A3 at (3 C decreases with increasing q and is very 
close to zero for q = 200. On a regular lattice the value is on the other hand closer 
to one, as the cluster boundaries are very complicated and the probability for two 
neighboring cluster to share a neighbor is very large. 

An interesting observation is that on dynamical triangulations both observ- 
ables seem to reach some kind of plateau value just before (3 C , then start to decrease 
again. This is more pronounced for A3. Looking at the observables for different 
lattice sizes indicate that their derivatives are discontinuous in (3 C . A plausible ex- 
planation of this is that these quantities approach constant values in the branched 
polymer phase, slightly larger than the limiting values W = 1 or A3 = 0. But as the 
model enters the pure gravity phase the values drop again, this time due to total 
magnetization. 

6 Discussion 

We have performed high statistics Monte Carlo simulations of the 10-state and 200- 
state Potts models coupled to quantum gravity. In a statistical mechanics notation 
this amounts to approximating the annealed average over triangulations of spherical 
topology. For both systems we observed a phase transition to a magnetized phase 
for sufficiently large (3, but the simulations suggest that the two systems belong to 
different universality classes. The transition for q = 200 is compatible with a third 
order transition, while data clearly favor a second order transition for q = 10. 

For the 200-state Potts model the dynamics at the transition agrees well 
the mean field calculations at q = 00. Just below the transition we observe many 
baby universes which are magnetized individually while the total magnetization is 
zero. Above the transition the spins of the individual baby universes are aligned. 
The data are not incompatible with the existence of a region II, below the critical 
point (3 C2 as suggested in fig.l. From the measurements of W, A3 and j s it seems 
that there may be a small region below (3 C2 where there are many magnetized baby 
universe and where j s ~ 1/2. This size of this region could well be comparable with 
the estimate (3 C2 — (3 Cl ~ 0.1 ± 0.05 mentioned in Sec. 2. 

Unfortunately we have not yet been able to locate the phase transition at 
(3 Cl . If there is a plateau where j s = 1/2 there has to be a second transition since it 
can been shown from general arguments that j s = l/(ra+ 1) if j s > 0, i.e. j s cannot 
change continuously to the value —1/2 for small (3. 

Maybe the most interesting result is that the critical behavior of the 10-state 
Potts model coupled to two-dimensional quantum gravity seems identical to that of 
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the 4-state Potts model. The 4-state Potts model has central charge c = 1, and the 
result indicates that there exists a whole family of g-state Potts models (probably 
up to q c > 10 of fig. 1) which on a regular lattice have first order transitions and no 
relation to central charge on a regular lattice, but which nevertheless after coupling 
to gravity lead to the same theory as a c = 1 conformal field theory and therefore 
should allow for a continuum interpretation, but in the case of the Potts models 
with 4 < q < q c only after coupling to gravity. Clearly this calls for a theoretical 
understanding, and it hints that the g-state Potts model coupled to two-dimensional 
quantum gravity might be exactly solvable. 
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